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We address the evolution of entanglement in bimodal continuous variable quantum systems interacting with 
two independent structured reservoirs. We derive an analytic expression for the entanglement of formation 
without performing the Markov and the secular approximations and study in details the entanglement dynamics 
for various types of structured reservoirs and for different reservoir temperatures, assuming the two modes 
initially excited in a twin-beam state. Our analytic solution allows us to identify three dynamical regimes 
characterized by different behaviors of the entanglement: the entanglement sudden death, the non-Markovian 
revival and the non-secular revival regimes. Remarkably, we find that, contrarily to the Markovian case, the 
short-time system-reservoir correlations in some cases destroy quickly the initial entanglement even at zero 
temperature. 

PACS numbers: 03.67.Mn,03.65.Yz 



I. INTRODUCTION 

Entanglement is an essential resource for quantum compu- 
tation and communication protocols [1]. However, this funda- 
mental quantum property is also fragile: the unavoidable in- 
teraction of quantum systems with their external environment 
leads to the irreversible loss of both quantum coherence (de- 
coherence) and quantum correlations in multi-partite systems 

A crucial requirement for a physical system to be of in- 
terest for quantum technologies is that the survival time of 
entanglement is longer than the time needed for information 
processing. Therefore it is important to develop a deep and 
precise understanding not only of the mechanisms leading to 
decoherence and entanglement losses but also of the dynami- 
cal features of these phenomena. Moreover, in view of recent 
developments in reservoir engineering techniques it is 

interesting to investigate situations in which decoherence and 
disentanglement can be controlled, for example through a pre- 
cise and accurate tuning of system and environment parame- 
ters. 

In this paper we consider the entanglement dynamics in 
noisy continuous variable (CV) quantum systems [7]. More 
specifically we focus our attention on a system of two non- 
interacting quantum harmonic oscillators bilinearly coupled to 
two independent structured reservoirs at temperature T. Our 
aim is to study the time evolution of the entanglement between 
the two oscillators for different temperature regimes, differ- 
ent system-reservoir parameters and different reservoir spec- 
tra. Rather than limiting ourself to present a plethora of dy- 
namical behaviors, we will try to identify general features in 
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order to single out universal properties of the disentanglement 
process, namely those properties that do not depend either on 
the specific model of reservoir chosen or on the specific ini- 
tial value of the entanglement. Moreover, we also compare 
the differences in the dynamics arising from different spectral 
distributions of the reservoir in order to identify those physical 
contexts leading to stronger or weaker entanglement losses. 

During the last decade numerous works dealing with losses 
and decoherence in bimodal CV quantum systems have ap- 
peared in the literature. In order to describe analytically the 
dynamics of such an open quantum system, approximations 
such as the Born-Markov and the secular (or rotating wave) 
approximations are typically performed [8, 9[ [T^[Tllll2l[T3ll . 
The Markovian approximation basically consists in neglect- 
ing the short-time correlations between system and reservoir 
arising because of the structure of the reservoir spectrum. 
This approximation is often performed together with the weak 
system-reservoir coupling assumption, also known as Born 
approximation. The Born and the Markov approximations are 
generally related. Indeed every time the coupling between the 
system and the environment is strong, and therefore the Born 
approximation is not appropriate, also the Markovian approx- 
imation cannot be consistently used. However, there exist 
situations of weak system-reservoir coupling and structured 
environment, where the system-reservoir correlations persist 
long enough to require a non-Markovian treatment, even in 
the weak coupling limit. In this paper we focus on these cases. 

We also investigate the validity of the secular approxima- 
tion (i.e. neglecting the counter-rotating terms in the Hamil- 
tonian) and we find that, even for weak couplings, a correct 
description of the short-time dynamics must take into account 
the nonsecular terms. In more detail, the validity of the secular 
approximation sensibly depends both on the reservoir temper- 
ature and on the system-reservoir parameters. 

Non-Markovian studies of bimodal CV quantum systems in 
a common reservoir have shown the existence of three differ- 
ent dynamical phases of the entanglement in the long time- 
limit, namely the sudden death, sudden death and revivals, 
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and no-sudden death phases ffl [l5ll . These phases depend 
not only on the system-reservoir parameters but also on the 
properties of the spectrum. In this paper we consider the case 
of two independent reservoirs and find a similar division in 
dynamical phases or regimes, namely, the entanglement sud- 
den death (ESD), the non-Markovian revival (NMRev) and 
the non-secular revival (NSRev) regimes. In our system, how- 
ever, the no-sudden death phase appears only at zero temper- 
ature and under very specific conditions. Moreover, we have 
discovered that the revivals may be due to two different phys- 
ical mechanisms, the non-Markovian finite reservoir memory 
or the presence of the non secular terms. 

Recent literature on non-Markovian CV dynamics, in the 
common reservoir scenario, includes Refs. Il6l [nil , while 
the independent reservoirs case was considered using a phe- 
nomenological approach in Ref. 01811 and using a numerical 
approach in Ref. 11191 l2(ill . In this paper, we extend in several 



directions the results we have obtained in Ref. |21|], where we 
limited our study to the high-T Ohmic reservoir in the secu- 
lar approximation. Here we solve the Master equation for our 
system without performing the secular approximation and in- 
vestigate quantitatively the entanglement dynamics using an 
analytic expression for the evolution of the entanglement of 
formation (EoF) 1I22I l23ll24ll . We assume the two oscillators 
initially excited in a twin-beam state (TWB, sometimes also 
referred to as two-mode squeezed states) and consider Ohmic, 
sub-Ohmic and super-Ohmic reservoirs at any temperature. 

The paper is organized as follows. In Sec. [n]we introduce 
the physical system, the Master equation and its general solu- 
tion through the characteristic function approach. In Sec. [Til] 
we review some preliminary concepts about two-mode Gaus- 
sian states and we present the general solution of the Mas- 
ter equation with an initial Gaussian state. We also introduce 
the TWB states, the concept of EoF for two-mode CV Gaus- 
sian states, and the types of reservoir spectra considered in 
the paper. In Sec. [[V]we present a detailed investigation on 
the validity of the secular approximation in our model. In 
Sec. |V] we discuss the dynamics of entanglement and an- 
alyze the three emerging dynamical regimes: ESD, NMRev 
and NSRev. Moreover, we give specific examples of the dy- 
namics of the EoF focusing on the high-T and low-T regimes. 
Finally, in Sec. [VT]we discuss and summarize our results, pre- 
senting conclusions and future prospectives. 



II. THE MASTER EQUATION 

We consider a system of two identical non-interacting quan- 
tum harmonic oscillators, each of them coupled to its own 
bosonic structured reservoir. The total Hamiltonian can be 
written as 

H = ^oaj% + X! E foikb] k b 5k 

3=1,2 j=l,2 k 

+ EE ^ + a \ ) ^ + b lk)> c 1 ) 

3=1,2 k 

with gjq the oscillators frequency, u)\k and L02k the frequencies 
of the reservoirs modes, a,j (a]) and bjk (b\ k ) the annihilation 



(creation) operators of the system and reservoirs harmonic os- 
cillators, respectively, and 7^ the coupling between the j-th 
oscillator and the fc-th mode of its environment. In the follow- 
ing we assume that the reservoirs have the same spectrum and 
are equally coupled to the oscillators. 

Since we are interested in the dynamics of the two oscil- 
lators only, we adopt a density matrix approach through the 
following local in time Master equation 1 25] 



■U(t)[X j ,[P j ,p(t)]} + -r(t)[X^p(t)}- 
i7(*)[^i»{^»P(*)}]. 



(2) 



where p(t) is the reduced density matrix, H° is the free 

Hamiltonian of the j-th oscillator, and Xj — (a,j + aj)/\/2 

and Pj = i(al — a,j)/y/2 are the quadrature operators. The 
effect of the interaction with the reservoirs is contained in the 
time-dependent coefficients of Eq. ©. The quantities A(£) 
and n(t) describe diffusion processes, j(t) is a damping term 
and r(t) renormalizes the free oscillator frequency ujq. 

It is worth noting that the Master equation (fJJ is exact, 
since neither the Born-Markov approximation nor the secu- 
lar approximation have been performed. The time dependent 
coefficients can be expressed as power series in the system- 
reservoir coupling constant. For weak couplings one can stop 
the expansion to second order and obtain analytic solutions 
for the coefficients. In the case of reservoirs in thermal equi- 
librium at temperature T, characterized by a spectral density 
J(oj), these expressions read 



A(t) 
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1] cos(ws) cos(wos), 



(3a) 



II(t) = cT / ds dojJ(uj)[2N(uj) + 1] cos(ws) sin(w s), 
Jo Jo 

(3b) 

rt r + oo 

7(t) = a 2 ds dujj(uj) sm(u>s) sin(uJos), (3c) 
Jo Jo 
pt r+00 

r(t) = a 2 ds dwJ(w) sin(o;s) cos(wos), (3d) 
Jo Jo 

where N(oj) — [exp(fiw/fc j BT) — is the mean number of 
photons with frequency u>, and a is the dimensionless system- 
reservoir coupling constant. 

By using the characteristic function approach 12611 . the so- 
lution of the Master equation © may be written as 



Xt (A) = e -^ T [W(t)®W(t)\K 



x X o{e- v{t)/2 [R-\t)®R-\t)]K), 



(4) 



w 3fc y u jk> 



where Xt(A) is the characteristic function at time t, xo is 
the characteristic function at the initial time t = 0, A = 
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(xt,px,X2,P2) is the two-dimensional phase space variables 
vector, r(i) = 2 f* j(t')dt', and W(t) and R-^t) are 2 x 2 
matrices. The former matrix is given by 



W(t) = e- r W[iJ- 1 (i)] T W(i)iJ~ 1 (<), 



(5) 



while the latter one, i?(t), contains rapidly oscillating terms. 
In the weak coupling limit R(t) takes the form 



R(t) = 



cos u>ot sin ojQt 
— sincjo^ cosa>o^ 



(6) 



Finally, W(t) = J * e r W M(s)ds with M(s) 
R T (s)M(s)R(s) and 

M w-l_n( s )/2 o 



The coefficient r(t) does not appear explicitly in the charac- 
teristic function solution because its contribution is negligible 
in the weak coupling regime [26]. The characteristic func- 
tion approach of Ref. 12611 is equivalent to other methods of 
solution of the Master equation (f2]i, as the Feynman- Vernon 
influence functional technique [27]. In this paper we use the 
former one because it allows to obtain an analytic solution in 
the weak coupling limit. 



If ci = C2 = the initial covariance matrix is block diagonal 
and the corresponding state is separable. Since each oscillator 
only interacts with its own environment, an initial separable 
state remains separable during all the evolution. For initial 
entangled states, however, the entanglement dynamics will in 
general depend on the initial value of the entanglement and 
on reservoir properties such as the spectral distribution, the 
temperature and the coupling constants. 

Since the evolution maintains the Gaussian character the 
evolved state is a two-mode Gaussian state with mean and co- 
variance matrix given by 



,-r(t)/2 



(R®R)%. 



( - 7 - ) Using Eqs. ©-([7]) we obtain 



Wt = 



- „-r(*) 



n( s ) 
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ds, 



where 



q u\ ( cos2dj t — sin2w t^ 



sin 2ojot — cos 2uJot I 



(11) 



a t = e~ r ^ (R © R)a (R © R) T + 2(W t © W t ), (12) 



(13) 



(14) 



III. ENTANGLEMENT DYNAMICS FOR GAUSSIAN 
STATES 

In this section we derive the explicit analytic solution for 
the characteristic function in the weak coupling limit already 
obtained in [26]. Remarkably, the evolution induced by the 
Master Equation (f2| corresponds to a Gaussian map, i.e. an 
initial Gaussian state maintains its character. It is thus possi- 
ble to obtain the expression of the covariance matrix at time t 
and then evaluate the EoF at any time for the two modes ini- 
tially excited in a TWB state. We also introduce the classes of 
spectral densities considered in the paper and show how the 
form of the time dependent coefficients ruling the dynamics 
can be simplified in the non-Markovian time scale. 



A. Analytic solution in the weak coupling limit 

Let us consider two-mode Gaussian states, i.e., those states 
characterized by a Gaussian characteristic function 



Xo(A) = exp j-iA T cr A - iA T X f 
We indicate with <tq the initial covariance matrix 



CT = 



A Co 
Cq B J ' 



(8) 



(9) 



where Ao = at, Bo = 61, Co = Diag(ci, cq), with a,b > 
and c%, cq real numbers, and 1 the 2x2 identity matrix. 
Moreover, 



Sa(«) 



sin 2uj$t cos 2uiQt 
cos 2wq£ — sin 2ojQt 



The covariance matrix at time t is given by 



At 






At 



with 

A t = A e~ T 



and 



Ct = 



A r + (A co -n si ) -(A si -n co ) 
-(A si -n co ) A r -(A co -n s 



ce r cos(2o;ot) ce r sin(2o;ot) 
ce~ r sin(2wot) — ce~ r cos(2o;oi) 



Xi n = Tr[p(0)(X 1 ,P 1 ,X 2 ,P 2 ) 5 



(10) 



A sl (t) 

n co (t) 
n si (t) 



- P -r(t) 



(15) 



(16) 



(17) 



(18) 



where we have introduced the function 

Ar(t) = e" r(t) f t e F ^A(s)ds (19) 
Jo 

and the secular coefficients 

A co (i) = e" r(t) / e r(s) A(s) cos[2w (t - s)}ds, (20a) 



f e r(s) A(s) sin[2cj - s)]ds, (20b) 

r(t) f e r ( s )n(s)cos[2w (<-s)]rfs, (20c) 
J o 

r(t) f e r( - s) U(s)sm[2uJo(t - s)]ds. (20d) 
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The explicit analytic expression of the coefficients above de- 
pends on both the reservoir spectral density and the tempera- 
ture. 

To further simplify the solution it is common to perform the 
so-called secular approximation. This approximation amounts 
at neglecting rapidly oscillating terms in the solution of the 
Master equation. In our case this means to assume that the 
coefficients (|20T > average out to zero. Stated another way, the 
secular solution is a coarse-grained expression of the exact 
one. In the next section we will critically examine the va- 
lidity of the secular approximation and derive the conditions 
of validity for different reservoir spectra and system-reservoir 
parameters. 

From now on we focus on TWB states, i.e., a set of 
Gaussian states whose covariance matrix (O has a = b = 
cosh(2r)/2 and c\ = — c 2 = sinh(2r)/2, with r > the 
squeezing parameter. Being pure states, their amount of en- 
tanglement is given by the entropy of entanglement Eo(r) = 
2 [cosh 2 r ln(coshr) — sinh 2 r ln(sinhr)] and, hence, it in- 
creases for increasing values of r. 



B. Entanglement of Formation 

A convenient and useful way of looking at the entanglement 
evolution in CV systems is by means of the EoF lf22l, 12311 . This 
quantity corresponds to the minimal amount of entanglement 
of any ensemble of pure bipartite states realizing the given 
state. In general it is not a simple task to derive an expression 
of the EoF for arbitrary states. Recently its expression for 
an arbitrary bimodal Gaussian state has been obtained in Ref. 
JH. 

We assume here that the initial state is a symmetric bipartite 
Gaussian state with covariance matrix given by Eq. (|9]). As we 
mentioned above, when this state interacts with two identical 
independent reservoirs, the Gaussian character is preserved 
and the evolved covariance matrix is given by Eq. ( fToT l. Due 
to the symmetry of the evolved state, the EoF is given by f23ll 



ln(x r( 



hi(x ri 



(21) 



withxm. = (k 2 _ + 1/4)/(2k_), k- = \J (a n - c+)(a n - c_) 
being the minimum symplectic eigenvalue of the CM er t , and 

(22) 




i? + ii-h±y/UT 



Ti-h) 



- {2hhf 



2Ji 

(23) 

where 1\ = dct[A t ], I3 = det [Ct] and I4 = dot [cr f ] are the 
symplectic invariants of <r t . Inserting Eqs. (fT~6b-(l20b into Eqs. 
(|2T1i-(|231 one obtains the analytic expression of the EoF for 
our system. 



C. Modeling the reservoir 

In order to obtain explicit expressions for the EoF we need 
to specify the properties of the bosonic reservoirs. We con- 



sider environments in thermal equilibrium at temperature T 
and we focus on the following class of Ohmic-like spectral 
distributions with an exponential cut-off function 



J s (w) 



1 'uj c 



(24) 



where u> c is the cut-off frequency. The case s = 1 corre- 
sponds to an Ohmic reservoir spectrum, characterized by a 
linear dependence on the frequency for ui -C u> c . For s > 1 
the spectrum is known as super-Ohmic while s < 1 describes 
a sub-Ohmic spectral distribution. For the sake of concrete- 
ness, in the following we consider the s — 3 super-Ohmic 
and the s = 1/2 sub-Ohmic cases. A more detailed discus- 
sion about the properties of these spectral distributions can be 
found, e.g., in Ref. JH. 

A closed form for the expressions of the time-dependent co- 
efficients given in Eq. (O can be obtained in the high-T and 
zero-T limits, i.e., for 2AT(w) + l m 2*gl and2N(w) + l w 1, 
respectively (See Appendix A). Therefore we focus on these 
two regimes. Inserting now the spectral distributions of Eq. 
(|24l . with s = 1, s = 1/2, and s = 3, into Eqs. © al- 
lows to determine the analytic form of the time-dependent 
coefficients. We notice that, in all three cases, after a time 
t w r c = w^ 1 , the coefficients attain their Markovian station- 
ary values and the system behaves according to the predictions 
of the Markovian theory. Here we are particularly interested 
in the non-Markovian short-time dynamics, and therefore we 
will focus on times t < t c . In this time interval, and in the 
weak coupling limit, we can expand the exponential terms ap- 
pearing in Eqs. $1% and (f20b in Taylor series. For example 
Eq. (T% becomes 



A r (t)~ / A(s)ds-T(t) f A(s)ds+ 
Jo Jo 

+ f T{s)A(s)ds + 0{a 4 ). 
Jo 



(25) 



Since A(t) oc a 2 and T(t) oc a 2 , in the weak coupling limit 
(a <C 1) and for short non-Markovian times the first term 
dominates and hence it is the only one that will be retained. 



IV. ON THE VALIDITY OF THE SECULAR 
APPROXIMATION 

In this section we question the validity of the secular ap- 
proximation by comparing the entanglement dynamics with 
or without the secular terms d20l ). As we will see, in general, 
the secular terms do influence the behavior of the entangle- 
ment in the short non-Markovian time scale. Depending on 
the value of certain parameters, however, the secular approx- 
imation in some cases turns out to give a good description of 
the dynamics. We have identified as main parameters influ- 
encing the time evolution the reservoir temperature, the pa- 
rameter x = uj c /ijjQ, and the form of the reservoir spectrum. 
We will discuss the effect of these parameters separately in 
the following three subsections dealing with the dynamics for 
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high-T reservoirs, T — reservoir, and with a comparison 
between different reservoir spectra. We will also consider if 
and how the validity of the secular approximation depends on 
the initial state of the system and, in particular, on the initial 
squeezing parameter r of TWBs. 

We recall that in this paper we will focus only on the dy- 
namics of entanglement. Different observables of the system 
may show different sensitivity to the secular approximation. 
Indeed we know that there exist a class of observables, e.g., 
the energy of the system oscillators, that are not influenced at 
all by this approximation 1 26] . 

A. High-temperatures regime 

We begin addressing the high-temperature limit ksT 3> 
Hluq, hio c , i.e., when the classical thermal energy fe^T is much 
larger than the typical energies exchanged in our system. In 
the following analysis we choose a temperature such that 
kBT/hu) c = 100, thus we can examine scenarios in which 
x = uj c /ujq > 0.1. 

We start analyzing the differences in the EoF evolution be- 
tween the secular result and the exact one in the case of an 
Ohmic distribution and x = 10. In Fig. [T]we plot the time evo- 
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creases and the secular coarse-grained solution describes well 
the dynamics of the entanglement in the short non-Markovian 
time-scale, independently from the initial condition. This be- 
havior is in agreement with the results of Ref. [28] where 
the weak coupling limit of the Master equation for quantum 
Brownian motion is discussed. In particular, in Ref. Il28ll it is 
shown that, in the high-temperature and weak coupling lim- 
its, the secular approximated Master equation is accurate only 
in the regime x <C 1, while for the other regimes the system 
behaves as if it were subjected to a squeezed reservoir. 




(b) 



FIG. 2: (Colors online) Comparison between the exact EoF dynamics 
(solid blue line) and the secular approximate dynamics (dashed red 
line) as a function of r = u> c t, with x = 0.2, (a) r = 1 and (b) 
r = 0.1. We set k B T/huj c = 100, a = 0.1. 

Summarizing, for high-T Ohmic reservoirs, the secular ap- 
proximation holds only in the regime x <C 1. This result is 
also valid for the sub-Ohmic and super-Ohmic environments. 



FIG. 1: (Colors online) Comparison between the exact EoF dynamics 
(solid blue line) and the secular approximate dynamics (dashed red 
line) as a function of t — uj c t, with x — 10, (a) r — 2 and (b) 
r = 0.5. We set k B T/huj c = 100, a = 0.1. 

lution of the EoF calculated using the secular approximated 
solution and using the exact solution in the regime i» 1 for 
two different initial TWB states. For both initial conditions 
the secular approximation fails. Remarkably, the exact solu- 
tion containing the nonsecular terms predicts a much longer 
disentanglement time. Furthermore the difference in the dis- 
entanglement time predicted by the exact and secular results 
increases for increasing values of initial entanglement, i.e., for 
larger values of r. This result is qualitatively independent of 
the analytic form of the reservoir spectrum, as we will see in 
Sec. II V CI were the effect of different spectra is considered. 

For intermediate values of the parameter x, x < 1, we ob- 
serve a stronger dependence on the initial value of entangle- 
ment. In Fig. [2] indeed, we see that for x = 0.2 and r = 0.1 
(small initial entanglement) the secular approximation works 
well, but for higher values of the initial entanglement, r = 1, 
nonsecular oscillations, absent in the secular approximated so- 
lution, are clearly visible in the exact dynamics. 

For x <C 1, finally, the nonsecular oscillations decrease in 
amplitude as the effective coupling with the environment de- 



li. Zero-temperature regime 

From previous studies on open quantum systems interacting 
with zero-temperature reservoirs we expect on the one hand a 
slower loss of entanglement |§] and on the other hand more 
pronounced non-Markovian features i30ll . with respect to the 
T ^ case. We will have a closer look at these general fea- 
tures of the dynamics in Sec. [V]and focus here on the validity 
of the secular approximation. 

We consider as an example a super-Ohmic reservoir with 
x = 0.3 and look at the dynamics of a TWB with a small 
amount of initial entanglement, r = 0.01. As shown in Fig. 
[3] the exact and the secular approximated dynamics sensi- 
bly agree in this situation. We have carefully examined the 
dynamical behavior for other values of x and of the initial 
squeezing parameter r reaching the conclusion that this is 
quite a general property of the system. Therefore, in the de- 
scription of bimodal CV quantum systems interacting with 
zero-T reservoirs, the secular approximation can always be 
preformed and the effect of the nonsecular terms is always 
negligible. This is a consequence of the fact that the secular 
terms (|20b are temperature dependent through the diffusion 
coefficients A(t) and U(t), and at T = their contribution is 
rather small. 
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FIG. 3: (Colors online) The exact (blue solid line) and the secular 
approximated dynamics (red dashed line) of the Ef vs r = u c t in 
a super-Ohmic reservoir at zero temperature for a — 0.1, r — 0.01 
and x = 0.3. The two curves almost overlap perfectly. 



C. Dependency on the reservoir spectrum 

To conclude our analysis of the secular approximation we 
look at the discrepancy between the secular and exact so- 
lutions for Ohmic, sub-Ohmic and super Ohmic reservoirs. 
Since in the zero-T case the secular approximation always 
works well, we focus on the high-T case and in particular on 
the x ^S> 1 regime, where the differences in the dynamics of 
the EoF are most pronounced. 

In Fig. [4] we compare the dynamics of Ep for the Ohmic, 
sub-Ohmic and super-Ohmic reservoirs as given by the sec- 
ular approximation, Fig. [4] (a), with the exact case, Fig. [4] 
(b), for x = 10 and r = 2. Comparing the two figures one 
clearly sees that the secular approximation does not affect the 
dynamics in an equal way for the three different spectral dis- 
tributions. The exact disentanglement time is almost tripled 
for the super-Ohmic environment and doubled for the Ohmic 
case. The sub-Ohmic case is less affected. In all cases, how- 
ever, and for all values of initial entanglement, the exact cal- 
culation predicts a longer survival time of entanglement. 




0.2 0.4 0.6 0.8 1.0 1.2 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 



FIG. 4: (Colors online) Comparison between the dynamics of Ef 
for the Ohmic (blue solid line), sub-Ohmic (red dashed line) and 
super-Ohmic (black dotted line) reservoir spectra using (a) the sec- 
ular approximated solution and (b) the exact solution in the high- 
temperature limit with ksT/huic = 100, a = 0.1, r = 2 and 
x = 10. 



V. ENTANGLEMENT DYNAMICS 
A. General features: Three different dynamical regimes 

In discrete variable quantum systems the phenomenon of 
ESD has recently received a lot of attention 03111 . In that con- 
text the basic system studied consists of two qubits interacting 
with either independent or common reservoirs. An exact solu- 
tion has been derived both for independent Ir32ll and for com- 
mon reservoirs 133J, and it has been shown that revivals of 
entanglement due to the reservoir memory may occur after an 
initial sudden death interval. The common reservoir scenario 
is generally characterized by a non-zero long time entangle- 
ment due to both the reservoir-mediated interaction between 
the qubits and the existence of a decoherence free subradiant 
state acting as an entanglement-trap ll34l HH . In the inde- 
pendent reservoir case, on the contrary, in the long time limit 
one always observes complete disentanglement. 

In the following we focus on the case in which the carriers 
of quantum information are continuous, instead of discrete, 
variable quantum systems. While the theoretical treatments 
of CV quantum channels is fundamentally different from the 
one of discrete channels, we find that some similarities in the 
entanglement dynamics do exist. In particular, in the com- 
mon reservoir scenario, due to the environment-mediated in- 
teraction between the two CV channels, the asymptotic long 
time entanglement maybe non-zero, even for high-T reser- 
voirs [ 14[ 1 1 311 . Moreover, non-Markovian studies show the 
occurrence of revivals of entanglement both in the common 
and in the independent reservoir cases HE! EH. 

Here we show that, for independent reservoirs and for 
x 3> 1, the phenomenon of ESD occurs both in the high- 
T and, for r -C 1, in the zero-T cases, independently from 
the reservoirs spectra, as one can see from Fig. Q] Fig. [4] 
and Fig. [5] (b). The ESD regime can directly be linked 
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FIG. 5: (Colors online) Dynamics of Ef at zero temperature for 
a = 0.1 in the case of Ohmic (blue line), sub-Ohmic (red line) and 
super-Ohmic (black line) reservoirs for x = 0.2 and r — 0.005 (a) 
and for x = 10 and r = 0.01 (b). 

to the behavior of the time-dependent coefficients appearing 
in the Master equation (|2). For high-T and x ^> 1, in- 
deed, independently from the reservoir spectra, the time de- 
pendent coefficients are always positive at every time instant 
H29I1 . In systems described by time-convolutionless Master 
equations non-Markovian features typically occur when the 
time-dependent coefficients temporarily attain negative values 
ll37i HH] . When this happens revivals of entanglement may oc- 
cur since the system restores partially the quantum coherence 
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previously lost due to the interaction with the environment. 
An example of non-Markovian revivals due to the reservoir 
memory effects, and therefore connected to the negativity of 
the time dependent coefficients, is shown in Fig. [6j where 
we plot the EoF dynamics for an Ohmic reservoir in the high- 
T limit. In general non-Markovian revivals of entanglement 
occur for x <C 1. In this case we have seen that the secular 
approximation works well and we know from previous studies 
(See Ref. IT29l,l39lo that the time dependent coefficients attain 
negative values for all reservoir spectra. 



E F 




FIG. 6: (Colors online) Ep dynamics vs r = u c t for an high-T 
Ohmic reservoir with fcsT '/huj c = 100, a = 0.1, x = 0.15 and 
r = 0.06. Note the presence of entanglement revivals due to negative 
values of the Master equation coefficients. 

Our exact approach allows us, moreover, to single out an- 
other type of revivals, not related directly to the negativity of 
the time-dependent coefficients, and therefore on the reservoir 
memory, but rather on the presence of nonsecular terms. An 
example of the nonsecular revivals is given in Fig. [7] where 
the exact and the secular solutions are compared. In this case 
the time-dependent coefficients are positive for t < 1.4 and 
the partial and temporary restoration of entanglement in that 
time interval is due to the presence of the counter-rotating 
terms in the microscopic Hamiltonian model. It is often be- 
lieved that the nonsecular or counter-rotating terms signifi- 
cantly affect the dynamics only in the strong coupling limit. 
For discrete variable systems indeed, a very recent study has 
shown the non-negligible effect of nonsecular terms in the 
strong coupling limit [40]. In this case the authors show that 
the exact dynamics causes a faster loss of entanglement with 
respect to the secular case. Here we show that also in the 
weak coupling limit these terms give a non-negligible con- 
tribution in the short non-Markovian time scale. In general, 
for intermediate values of x, the dynamics will display both 
non-secular and non-Markovian revivals. Indeed, if we look 
at the dynamics for r > 1.4 we see that a non-Markovian 
revival, with a superimposed nonsecular revival, occurs at 
1.5 < t < 2.5, as shown in Fig. [7] (b). The border be- 
tween the NMRev and the NSRev dynamical regimes is there- 
fore a blurred region in which both effects occur at the same 
time. In this case both the reservoir memory and the nonsec- 
ular terms contribute to the reappearance of previously lost 
entanglement. 




T 



T 



FIG. 7: (Colors online) Ef for a Sub-Ohmic reservoir in the high-T 
limit (k B T/hio c = 100) with a = 0.1, r = 2 and x = 0.3. The 
solid blue line is the exact solution while the dashed red line is the 
secular solution. The inset is a magnification of the temporal region 
< r < 1.5. 



A more detailed discussion is required to explain the pres- 
ence of ESD in the zero-temperature cases shown in [5] (a)- 
(b). The Markovian theory of two-mode continuous variable 
channels predicts, for both the common and the independent 
reservoirs, the existence of a finite time of disentanglement 
for an initial TWB state when T > 18[] . In the independent 
reservoirs model at T — 0, however, the Markovian disen- 
tanglement time is infinite (no ESD). Since the exact Mas- 
ter equation (f2| coincides with the approximate Born-Markov 
Master equation for weak couplings and for times longer than 
the reservoirs correlation time, one would expect our non- 
Markovian model to give the same prediction for the disen- 
tanglement time than the Markovian one. Stated another way, 
one would not expect ESD. 

However, one should keep in mind that the Markovian ap- 
proximation is always a coarse graining in time and therefore 
it does not allow us to predict the short time non-Markovian 
behavior. If at short times the initial entanglement is lost and 
no non-Markovian revivals occur, entanglement cannot reap- 
pear at longer times. Consequently the entanglement will re- 
main zero also in the asymptotic Markovian long-time region. 
This is exactly what may happen when the initial amount of 
entanglement is small (r <C 1). In this case, indeed, for 
some reservoirs spectra and values of x, the short time non- 
Markovian dynamics shows the occurrence of sudden death 
of entanglement. Since the state remains separable for times 
greater than the reservoirs correlation time, entanglement re- 
vivals cannot appear. 

For higher values of initial entanglement, on the other hand, 
the exact non-Markovian theory does not lead to a sudden 
death in the short time scale, thus the state is still entangled 
when reaching the Markovian time-region and therefore the 
Markovian prediction of an infinite disentanglement time at 
T = still holds. 

Summarizing, for x 3> 1 ESD occurs, independently from 
the reservoir spectrum, both in high T reservoirs (for all val- 
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ues of r) and in zero-T reservoirs (for r <C 1). When x <C 1 
we are generally in the NMRev region, independently from 
the reservoir spectrum. One should note, however, that if the 
initial entanglement is very small (r -C 0.1) entanglement 
oscillations do not have time to take place and only ESD is 
observed. For intermediate values of x the dynamical regimes 
strongly depend both on the reservoir spectrum and on the 
initial entanglement. More specifically, for high-T reservoirs 
one can have any of the three ESD, NMRev and NSRev be- 
haviors, as well as a combination of NMRev and NSRev. For 
zero-T reservoirs the ESD or NMRev regimes exist only if the 
initial amount of entanglement is small (r <C 1), for other ini- 
tial values of r entanglement is never lost, in accordance with 
the Markovian theory. 



B. Comparative study of Ohmic, sub-Ohmic and super-Ohmic 
reservoirs 

In this section we investigate the differences in the loss 
of entanglement due to different reservoir spectra. Different 
physical systems are characterized by different environmental 
spectral densities, e.g., it is well known that solid-state sys- 
tems are subjected to sub-Ohmic 1// noise. Such a compara- 
tive study, hence, allows to understand which physical context 
is more 'quantum information friendly', in the sense of allow- 
ing entanglement to live longer. 

We begin considering the high-T reservoir case. In Fig. [4] 
(b) we have seen that in the ESD regime and for high tem- 
peratures the behavior of entanglement is qualitatively simi- 
lar. The disentanglement time is not strongly dependent on 
the reservoir spectrum and the sub-Ohmic environment dis- 
plays the faster loss of entanglement. For intermediate values 
of x, however, the time evolution of Ep shows a much richer 
behavior and a much stronger dependence on the form of the 
spectrum, as one can see in Fig. [8] In this case the super- 
Ohmic environment shows a much faster loss of entanglement 
than the Ohmic and sub-Ohmic. In the Ohmic and sub-Ohmic 
cases one can clearly see the nonsecular oscillations superim- 
posed to the non-Markovian oscillations in the dynamics of 
Ep, the latter ones having longer period and larger amplitude. 

The entanglement dynamics for zero-T reservoirs is plot- 
ted in Fig. [5] for (a) x = 0.2 and (b) x = 10. Also in this 
case the super-Ohmic reservoir shows a much faster entangle- 
ment loss than the sub-Ohmic and the Ohmic ones. This is 
especially evident in the intermediate x region of Fig. 0(a). 
The Markovian and RW theory of entanglement dynamics for 
TWBs in independent reservoirs predicts that the disentangle- 
ment time (or separability time) should approach infinity for 
T — > iU]. Our results demonstrate that this conclusion, in 
the case of small initial entanglement, is a consequence of the 
Markovian approximation and that the exact non-Markovian 
theory predicts that, even for weak system-reservoir coupling, 
in the x ^S> 1 region the entanglement survives only for a 
short time. For intermediate values of x, however, the dis- 
entanglement time approaches the Markovian prediction for 
the sub-Ohmic and Ohmic cases. Indeed we see that after ini- 
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FIG. 8: (Colors online) Dynamics of the Ef in the high-T temper- 
ature limit kBT/luu c = 100 for r = 2 and x — 0.2 in the case 
of Ohmic (blue solid line), sub-Ohmic (red dashed line) and super- 
Ohmic (black dotted line) environments. 

tial non-Markovian oscillations Ep approaches its stationary 
non-zero Markovian value. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have studied a bimodal CV quantum sys- 
tem interacting with independent structured reservoirs in ther- 
mal equilibrium. We focus on the dynamics of the entan- 
glement, as measured by the entanglement of formation, for 
the two modes initially excited in a twin-beam state and for 
different reservoir spectral distributions (Ohmic, sub-Ohmic 
and super-Ohmic). Under the only assumption of weak cou- 
pling we have obtained an exact analytic solution for the time- 
dependent two-mode covariance matrix describing the state of 
our system in the short time non-Markovian limit. 

In the first part of the paper we unravel the role of the sec- 
ular approximation in our specific system and in particular in 
the context of the entanglement dynamics. By comparing the 
exact solution with the solution in the secular assumption we 
found that, in general, the high-T dynamics is affected by this 
approximation, while the T = case is not. More specif- 
ically at high temperatures we have shown that for x 3> 1 
both the exact and secular solutions predict the occurrence of 
ESD. However the entanglement persists longer in the exact 
solution. On the contrary in the limit of x <C 1 the solutions 
agree. These two results are independent from the analytic ex- 
pression of the reservoirs spectral distributions and from the 
initial amount of the entanglement. For intermediate values 
of x the situation is more complicated, and the validity of the 
secular approximation strongly depends on the expression of 
the spectral distribution and on the initial amount of entangle- 
ment. 

The second aim of the paper was the investigation of the 
entanglement dynamics as a function of the reservoir spec- 
trum, the temperature and the initial amount of entangle- 
ment. Essentially we observed the presence of three differ- 
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ent behaviors: sudden death of entanglement (ESD), non- 
Markovian entanglement revivals and oscillations (NMRev) 
and oscillations or revivals related only to the secular coeffi- 
cients (NSRev). 

At high temperatures ESD appears for each value of the ini- 
tial entanglement. Moreover, for x ^S> 1 there are no revivals 
while they are typical of the dynamics for x <C 1. Because 
in this limit the secular and exact dynamics almost coincide, 
these revivals are due to the negativity of the Master equa- 
tion coefficients (NMRev). For intermediate values of x, the 
time evolution strongly depends on the initial amount of en- 
tanglement. If the entanglement is small (r < 1) only ESD 
is observed. For larger value of r, in general, the behavior is 
characterized by oscillations and revivals. Some of these re- 
vivals occur in correspondence of positive value of the time- 
dependent coefficients. Therefore they exist as a consequence 
of the secular coefficients only (NSRev). 

The situation for T = is characterized by a slower rate 
of entanglement deterioration. Therefore, for short times, en- 
tanglement sudden death and revivals can be observed only 
for very small initial entanglement (r <C 1). In these cases, 
and when x 3> 1, ESD exists independently from the reser- 
voir spectrum while for x <C 1 EDS is present in the super- 
Ohmic case only. Hence, the asymptotic long time Markovian 
dynamics of entanglement, and therefore also the Markovian 
prediction about the disentanglement time, may be strongly 
affected by the non-Markovian short time correlations. When 



this happens, the non-Markovian theory predicts a finite dis- 
entanglement time in contrast to the Markovian prediction. 
When r > 0.1 the short time non-Markovian dynamics is 
characterized by oscillations only, the EoF remaining posi- 
tive. Therefore, for long times the Markovian prediction of an 
infinite disentanglement time is recovered. 

In recent years there have been a lot of interest in the entan- 
glement dynamics in C V quantum channels, both for common 
and independent reservoirs. Our work finds its place in this 
context as an attempt to investigate the non-Markovian short 
time dynamics of entanglement in different physical scenar- 
ios. We believe that our results, showing the effects of dif- 
ferent reservoirs on the time evolution of entanglement in CV 
quantum channels, will pave the way to the implementation 
of engineered reservoir control schemes as the one recently 
reported in Ref. 14111 for qubits. 
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APPENDIX: TIME DEPENDENT COEFFICIENTS AT THE SECOND ORDER IN a 



Here we provide the exact analytic expressions of the time-dependent coefficients of the Master Equation given in (01. In 
each subsection we consider a single reservoir spectral function and evaluate the temperature independent damping coefficient 
j(t), the diffusion coefficients in the high temperature regime Ar(i) and ITr(i), and the diffusion coefficients at T = 0, Ao(i) 
and Ho(t). The expression for T(t), Ar(i) and the secular terms ( 1201 follow through. We made use of the following special 
mathematical functions 1,4211 
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3. Super-Ohmic Reservoir s = 3 
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